Singular perturbation near mode-coupling transition 



Mami Iwata and Shin-ichi Sasa 

Department of Pure and Applied Sciences, University of Tokyo, Komaba, Tokyo 
153-8902, Japan 

Abstract. We study the simplest mode-coupling equation which describes the time 
correlation function of the spherical p-spin glass model. We formulate a systematic 
perturbation theory near the mode-coupling transition point by introducing multiple 
time scales. In this formulation, the invariance with respect to the dilatation of time in 
a late stage yields an arbitrary constant in a leading order expression of the solution. 
The value of this constant is determined by a solvability condition associated with a 
linear singular equation for perturbative corrections in the late stage. The solution 
thus constructed provides exactly the a-relaxation time. 



1. Introduction 

About a quarter of a century ago, a peculiar type of slow relaxation was discovered 
theoretically in researches of glassy systems [H [2]. This relaxation behavior is 
characterized by two different time scales, both of which diverge at a temperature. 
An illustrative example exhibiting such a behavior is the spherical p-spin glass model 
[3]. The normalized time-correlation function of total magnetization in this model 
turned out to satisfy exactly the so-called mode-coupling equation, which is written as 



for the case p = 3. Here, the initial condition is given by 0(0) = 1, and the parameter g 
is proportional to the square of the inverse temperature. Since this equation is derived 
under an assumption that the system possesses the stationarity, is valid only in a 
regime < g < Qc, where gc will be given later. 

A remarkable feature of ([T]) is the existence of nonlinear memory. One can regard 
([1]) as one of the simplest equations that characterize a universality class consisting 
of models with nonlinear memory. Indeed, some qualitatively new features in glassy 
systems have been uncovered by studying ([1]) and its extended forms. (See Ref. [1] as 
a review.) In particular, two divergent time scales were found just below the transition 
point g = gc, and the precise values of the exponents characterizing the divergences were 
determined. Furthermore, extensive studies have been attempted so as to construct the 
solution in a systematic manner. On the basis of past achievements, in the present 
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paper, we propose a perturbation method for analyzing ([T]), which might shed new hght 
on the nature near the mode couphng transition. 

We shall address the question we study in this paper. Let foo be the value of 
(j)(t oo). We substitute = G{t) + f^o into ([1]) and take the limit t — oo. We then 
obtain 

-/oo + ^?/^(l-/oo)=0, (2) 

where we have used the relation G(0) = 1 — f^o- From the graph gf^(l — foe) as a 
function of foo, we find that the non-trivial solution (/oo 7^ 0) appears when g > gc = 4:. 
This transition is called the mode-coupling transition. Note that foo = 1/2 when g = g^. 
Below we express this value of foo as fc- We then introduce a small positive parameter 
e by setting g = gc — e, and we denote the solution of by In this paper, we 

formulate a perturbation theory for ([1]). As a result, we obtain an asymptotic form of 
(j)e{t) in the small e limit. 

More concretely, for a given small positive e, we want to express (f)e{t) in terms of 
e and e-independent functions. For readers' reference, in figure [H (left), we display the 
numerical solution (/>e(t) with e = 10~^. Here, when solving ([1]), we used the algorithm 
proposed in Ref. [5]. In figure [1] (right), we also display the e-dependence of the a- 
relaxation time Tq, which is defined by (pei^a) = 1/4. We want to calculate Tq, based on 
our theory. 




Figure 1. (j>e{t) with e — 10"'^ (left), a-relaxation time as a function of e (right). 
The circle symbols represent the result of numerical simulation of ([T]) . The dotted line 
corresponds to the theoretical calculation Tq = 20e~^'^^ given by (|43| . 



This paper is organized as follows. In section [21 we set up our theory. In particular, 
we give a useful expression of a perturbative solution. This section also includes a review 
of known facts in order to have a self-contained description. In section [3], we formulate 
a systematic perturbation theory on the basis of our expression, and we determine a 
leading order form of the solution. We check the validity of our theory by comparing 
our theoretical result of Tq, with that measured by direct numerical simulations of 
([1]). Section m is devoted to remarks on possible future studies. Technical details are 
summarized in Appendices. 
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2.1. Solution with e = 



We first investigate the solution (poit). It is expressed by (f>o{t) = Go(t) + fc with a 
function Go(^) which decays to as t ^ oo. The equation for Go is written as 



dtGo{t) + f, + Goit)+g, 



fds{f, + Go{t-s)fdsGo{s) = 0. 
Jo 



(3) 



An asymptotic form of Go{t) in the large t limit can be derived by employing a formula 
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for any x < 1 and t > 0, where T{x) is the Gamma function. The result is 
Go(t) ~ cot"^ 

where a is a constant that satisfies a relation 
r^fl-a) 1 
2' 



(4) 



(5) 



(6) 



r(l -2a) 

The value of a is estimated as a = 0.395. In figure [2], we display the graphs of 4>o{t) and 
Go(^)) which are calculated numerically. 



^app 



As shown in Appendix A, an approximate expression of Cq is calculated as Cq — 
(a/(l — a))°'{l — a)/2°'~^'^ by a matching procedure. Its value, 0.194, is not far from 
Co = 0.25 obtained from a numerical fitting of the graph of Go{t). It might be possible 
to improve the approximation in a systematic manner. However, in this paper, we do 
not pursue such improvements. The important thing here is that the e-independent 
function Go(^) is defined with understanding of its asymptotic form. 
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Figure 2. (j)o{t) (left) and Go{t) in the log-log plot (right). The dotted line in the 
right figure represents 0.25i~°'^^^. 



2.2. Expression of the solution with e > 

One may expect that the solution 0e(^) is close to 0o(^)- However, recall that (poit — > oo) 
changes discontinuously from to 1/2 when g passes at g = gc from below. This 
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fact means that a small perturbation from g = g^ {e = 0) yields a singular behavior. 
Examples of such singular perturbation can be seen in Refs. [6, 7J. 

Before formulating effects of the perturbation, we conjecture a functional form of 
the solution First, should be close to 4>o(t) in an early stage where 0^ > fc- 

Since 4>o{t) — > /c as t ^ cxd, the trajectory (peit) stays in a region near fc for a long 
time. However, since there is no non-trivial solution /oo 7^ for positive e, 0^ goes away 
from the region 0^ ~ and finally approaches the origin 0^ = 0. Such a behavior is 
substantially different from 0o(t). We thus introduce a quantity A{e'^'H) that describes 
the relaxation behavior from 0^ ~ to 0^ = 0, where ^(0) = fc- The functional form 
of A is independent of e, while its argument is the scaled time t2 = e'^^t with a positive 
constant 72. We also expect that a switching from 0e(t) ~ 0o(^) in the early stage to 
0e(t) ~ A{e^^t) in the late stage occurs around another characteristic time of 0{e~^^) 
with a positive constant 71. Keeping this behavior in mind, we express the solution as 



where the switching function satisfies 0(0) = 1 and 0(oo) = 0. The functional form 
of is independent of e, while its argument depends on e. y:>e(t) represents a small 
correction that satisfies (/^^ — in the limit e — > for any t. 0, A, 71, 72 and ife{t) will 
be determined later. 

In order to have a simple description, we define a set of scaled coordinates (to, ^1,^2) 
on the time axis as ti = e'^H, where 70 = 0. Throughout the paper, a time coordinate 
with an integer subscript represents the scaled coordinate determined by the subscript. 
Note that to, ti and t2 appear as the arguments of Gq, © and A, respectively. Physically, 
the first relaxation occurs in the early stage to ~ 0{1)] the behavior around 0e = fc is 
observed in the intermediate stage ti ~ 0(1); and the relaxation behavior from 0^ ~ fc 
to 0e = is described in the late stage t2 ~ 0(1). In researches of glassy systems, the 
intermediate and the late stages are termed the /^-relaxation regime and the a-relaxation 
regime, respectively. 

2.3. Equation for A 

We substitute ([7]) into ([I]) and take the limit e — with the scaled time ^2 = e^^t fixed. 
We then obtain 



In this paper, the prime symbol represents the differentiation with respect to the 
argument of the function. The equation ([8]) provides the explicit definition of A with 
the condition A{0) = fc- However, this equation cannot determine A uniquely. Indeed, 
for a given solution ^(^2) of ([8]), A{Xt2) with any positive A is another solution of ([8]). 
This dilatational symmetry is a remarkable property of ([8]). 

Here, by analyzing the short time behavior in ([8]), one can confirm that ^(^2) — fc 
is proportional to t2 in the small ^2 limit. Thus, we can choose a special solution of 
A such that v4'(0) = —1. In the argument below, A represents this special solution; 



0,(t) = Go(t)0(e^^t) + A{e^H) + y;,(t). 



(7) 
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and the other solutions are described by A{Xt2). For later convenience, we define Ax by 
A\{t2) = A{Xt2), and A in ([7j) is replaced with Ax- In particular, we have 

Axit2) = fc- Xt2 + 0(t2) (9) 

in the limit t2 ^ 0. Note that A is an arbitrary constant until a special requirement is 
imposed. The functional form of A{t2) can be obtained by solving ([8]) numerically with 
a simple discretization of time. We display the graph of ^(^2) in figure [3l It should be 
noted here that the mathematical determination of the functional form is not the heart 
of the problem. The important thing is that the e-independent function A is defined 
without any ambiguities. 




t t 



Figure 3. A{t) (left) and its semi- log plot (right). 



We express the dilatational symmetry in terms of a mathematical equality. Let us 
define 

$2(^2; Ax) = Ax{t2) - gdl - fc)Al{t2) +gc r ds2Al{t2 - s2)A'^{s2).m 

Jo 

Since ([H]) is identical to $2(^2;^) = 0, the dilatational symmetry is expressed by 
$(^2; ^a) = for any A. Then, taking the derivative with respect to A, we obtain 

dS2LA{t2,S2)dxAx{s2) =0, (11) 

where is the linearized operator around Ax, which is defined by 

L4t.s,)^'-M2^. (12) 

(>Ax{S2) 

Its explicit form is given by 

LA{t2, S2) = 6{t2 - S2) (1 - 2g,{l - f,)Ax{s2) + gJl) 

+ gc0{t2 - S2)2A'x{t2 - S2) {Ax{s2) + Ax{t2 - S2)) . (13) 

It will be found below that (ITTl) plays a key role in our formulation. 

Now, from the definition of 71, we have /c + G'o(e~'^^) ^ ^^(e'^^"'''^). By substituting 
([5]) and ([9]) into this relation, we obtain 

71 = T^. (14) 
1 + a 
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2.4- Functional form of & 

We can formulate a systematic perturbation theory with employing an arbitrary 
switching function 0(ti) when it decays faster than a power-law function as 
we will see in the next section. For example, one can choose a physically reasonable 
form 



where t = t^.e""'^ corresponds to the time when the graph 0o(^) is closest to that of 
Ax{e"''^t). That is, t^ satisfies acot~"~^ = A. Note, however, that there is no reason 
that we must choose this form. Indeed, other forms such as 0(ti) = exp(— 2ti/i(:c) and 
0(ti) = exp(— (ti/tc)^) might also be physically reasonable to the same extent as (fT5|) . 
Of course, the final result should be independent of the choice of the functional form. 

2.5. Summary 

In our formulation, we set the unperturbative solution 0u as 



Go and A are already determined. 71 is connected with 72 in (1141) . 6 is assumed to 
take an arbitrary form. Thus, the problem we solve is the determination of 72 and A as 
well as the perturbative calculation of the correction ^Peif)- Note that 72 and A appear 
in the leading order expression of the solution. In particular, since the value of A has 
never been known, the calculation of A is a cornerstone of our theory. 

3. Systematic perturbation 

3.1. preliminary 

For any trajectory ip{t) with V'(O) = 1, we define 



e(ti) = exp(-ti/tc), 



(15) 



0u(t) =G'o(t)e(e^^t)+A,(e^^t), 
and we express the perturbative solution 0^ by 

0,(t) =0,(t)+(^e(t). 



(16) 



(17) 




(18) 



Let L^{t, s; ip) be the linearized operator of F^{t] ip), which is defined by 




(19) 



Its explicit form is written as 

L,{t,s;tlj)=6{t-s){l + g) + 5'{t-s) 

- ge{t - s) {dsi^\t - s) - 2i;{sW{t - s)) . 
The mode-coupling equation ([1]) is expressed by 

Feit- 0,) = 0. 



(20) 



(21) 
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3.2. Calculation 

The substitution of ( fTTl) into ( pTl) yields non-trivial e dependences through the evaluation 
of the integration term in ( |T8i) with using the scaled coordinates. In order to avoid 
a complicated description, we focus our presentation on an important part of the 
calculation. 

We first evaluate F^{e~"'^t2:,(t)u) in the small e limit with t2 fixed. As explained in 



Appendix B we derive 

F^-^H^, 0u) ^ e-A(t2) + e^^-^^(i-'^)2(?e(/c + A^{t2))A'^{t2)c,e, (22) 
9c 

where higher order terms of 0{e''^~'^^^^~'^°''>) are neglected, and 6' is a constant determined 
by 

/■oo 

e= dss-^eis). (23) 

Jo 

We here make three important remarks on ( l22l) . First, if we did not introduce the 
switching function in the expression of the solution ([7]), we would have a form rather 
different from (!22|) . for which the analysis seems to be hard. Second, the function G 
should provide a finite value of 6. This means that 0(ti) should decay faster than 
a power-law function Third, the two terms in fl22l) should balance each other. 

Otherwise, a contradiction occurs. (See an argument below ( 1371) .) The last remark leads 
to the relation 72 — 71(1 — 0) = 1- By combining it with (fT4l) . we obtain well-known 
results 

and 

+ i (25) 

which correspond to the exponents characterizing divergences of the /^-relaxation time 
and the a-relaxation time in glassy systems, respectively. Then, fl22|) with flMl) and fl25l) 
becomes 

F,{e-^H2; 0u) = e4'\t2) + 0{e'/'), (26) 
where JFg^"* is the e-independent function given by 

4'\t2) = -Ax{t2) + 2g,{f, + Ax{t2))A'^{t2)coe. (27) 
9c 

Furthermore, we can prove 

e-^'L,{e-^H2, e'^^Ss; 0u) = LA{t2, S2) + 0(e), (28) 

where La is given by f|T3|) . 

Now, let us compare (l26l) and ( 12T]) with f |T71) . We then find that the perturbative 
correction is expressed as 

V>e{e-^H2)=e^S\t2) + Oie'/') (29) 
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in the regime ^2 ~ ^(1) with the hmit e ^ 0. We also write 

^,(e-^Hi)=e"v5l")(ti) + o(e") (30) 
in the regime ti ~ 0(1) with the hmit e ^ 0, where a is a positive constant. Then, 



(121]) with ([ilD becomes 

POO 

Feie-^H^; K + e"<^i"^) + e / ds2LA{t2, S2)^^2\s2) = o(e), (31) 

where the contribution of v5e(to) is included in the right-hand side. Here, by an argument 
similar to Appendix B , we can estimate 

F,(e-^^t2; 0u + e"^i"^) = eT!^\t2) + 0{e^^"^). (32) 

In order to describe a theoretical framework in a simple manner, for the moment, we 
focus on the case a > 1/2. The other case a < 1/2 will be discussed in section |331 The 
equation for (^2^''(t2) is then simply written as 

POO 

/ dS2LA{t2,S2)^^2\s2) = -:F^2\t2). (33) 

Jo 

3.3. Solvability condition 

We notice that (!33|) is a linear equation for (^2^\ which is singular because there exists 
the zero eigenvector $o = d\A\ associated with the dilatational symmetry (ITTl) . Let 
be the adjoint zero eigenvector that satisfies 



ds2LA{s2M)%{s2) = Q. (34) 
Then, there exists a solution of ( |33l) only when the condition 

dt2^l{t2)J'^^\t2) = (35) 



is satisfied. Otherwise, (133|1 leads to 7^ and hence there is no solution ^p^2^ ■ The 
equality ( !35l) is called the solvability condition for the singular equation (!33ll . Note, 
however, that the solvability condition is not satisfied as an identity. Here, let us recall 
that A is still an arbitrary constant. Thus, we are allowed to determine the value of A 
so that the solvability condition ( l35l) can be satisfied. Only for this special value of A, 
the perturbation theory can be formulated consistently. 

Concretely, since we find that $0(^2) = 5(s2) from ( !34l) with (fT3l) . the solvability 
condition fl35l) becomes 

J^^^^(0)=0. (36) 

The explicit form JF2^''(0) obtained from (127|) leads to 



Here, let us go back to fl22]) . If 72 — 71 (1 — a) were not equal to 1, the condition 
could not be satisfied for any positive A. In this sense, one can regard that the solvability 
condition determines the exponent 72 as well as the constant A. 
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3.4- Determination of \ 

Apparently, ( l37l) shows that A depends on the choice of 6. However, the situation is 
a httle bit comphcated. We shall explain the way how to determine the value of A in 
detail. 

We study the case e ^ with ti fixed. In this limit, can be expressed as 

Ue-^H,) = /e + e^/2[cQt--e(ti) - Ati] + ^.{e-^H,), (38) 

where we have used ([5]) and ([9]). Since 6 is arbitrary, {p^{e~'^'^ti) includes a term of 
0(e^/^) unless a special 9 is employed. This means that a in ( 130|) is equal to 1/2. As is 
seen from (1321) . when a = 1/2, Lp^\t2) must be calculated with taking account of ^^i^'^^- 
Concretely, the right hand side of (1331) should be replaced with —J^2^\t2), where 

F.{e-^H2; 0. + e^/^^S'^'^) = + o(e). (39) 

Then, the solvability condition (|36|) is also replaced with ^2^\o) = 0. Without the 
replacement, ( 1371) provides nothing more than an approximation of A. For example, (1371) 
with (US]) leads to A = [(64r(l-a))(i+'^)/(2a)^i/«^(i-a)/(2a)]-i _ q as one approximate 

value. 

Now, let us calculate the precise value of A. One natural method is to choose a 
functional form of G so that the condition a > 1/2 is satisfied. We denote this special 
by B*. Then, for a given B, the correction <^[^'^^^ is determined by 

cotr"0(^i) + = cot^^e.iti). (40) 

Therefore, (j36l) using B* is equivalent to ^2^\o) = using B. In other words, the 
precise calculation of A starting from B can be done through B*. This also indicates 
explicitly that the final and precise result is independent of the choice of B. In any 
cases, the problem is focused on the calculation of B*. 

As explained in [Appendix C[ we can derive the equation for Q(ti) = Coti"'Q^:(ti) in 
the form 

^-8A r dsi[Q{s,)-Q{t,)/2] 
o Jo 



+2Q^ + 4 / dsilQit, - s,) - Q(ti)]g'(si) = 0. 
Jo 



(41) 



We study this equation by regarding A as a parameter whose value is not specified 
beforehand. We denote this solution by Q{ti;\). For almost all A, Q(ti; X) is not 
bounded as ti — oo, while there exists the special value A,,, such that Q(ti; A*) — * as 
ti ^ oo. A necessary condition for this property is easily derived by considering the 
limit ti ^ cx) in ( HTI) : 

1 -1 



Aj, — — 
64 



(42) 



This is equivalent to the expression fl37l) that determines the value of A by the solvability 
condition fl5^ under the assumption a > 1/2. Therefore, once we find A^, such that 
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Q(ti;X*) — s> as ti oo, this A* is the precise value of A that we want to have. 
Simultaneously, we obtain 0*(ti) from A*). 

The problem of finding A* is investigated by a shooting method. We first solve 
(14T|) numerically for a given A. Basically, we employ a simple discretization method. In 
order to treat properly the singular behavior near t = 0, we utilize the result of the short 
time expansion of Q{ti; A) near t = 0. (See Appendix D for the short time expansion.) 
Suppose that we already investigated the system with A^, k = 0,1, ■ ■ ■ ,n. We here note 
that Q{ti; A) — oo as ti ^ oo when A = and that Q{ti; A) ^ oo as ti ^ oo when 
A is sufficiently large. Based on this observation, we define /i^ = maxA^ such that 
Q{ti] Afe) — oo as ti oo, and /!„ = min Afc such that Q(ti; Xk) ^ oo as ti ^ oo. 
We then choose A„+i as A„+i = (/i^ + ^1^)12. Starting from Aq = and Ai = 1, we can 
determine the sequence {A^} for which Aoo = lini„^oo A„ exists. From the construction 
method, Q{ti; Xoo) — as ti — > oo. Therefore, A^, is given by A^o- By performing this 
procedure numerically, we estimate A,. = 0.017. In this manner, we have determined the 
precise value of A and the function 0^. We display the functional form of 0^, in figure 

SI 



0* 




0* 






100 



Figure 4. 6,(t) (left) and its semi-log plot (right). 



3.5. Remarks 

At the end of this section, we make two remarks. First, as a demonstration of our result, 
we study the a-relaxation time Tq defined by (peija) = 1/4. Let ta be A{ta) = 1/4. 
Then, from the expression of the solution (fT7|) . is estimated as = (r^/A)e~'^2. By 
using the value ta = 0.346 obtained from numerical integration of (IE]), we arrive at the 
theoretical prediction 

= 20e'^-'^\ (43) 

In figure [1] (right), we display the result of numerical simulations of ([T]) with e = 0.1 x2~-', 
j = 0, ■ ■ ■, 10. The numerical data are perfectly placed on the theoretical result (H3|) . 
This is an evidence that the expression (1371) is correct. 

The second remark is on the systematic formulation. In principle, higher order 
terms such as (p f^\t,) and <^5-^^H^.) 

can also be calculated in a manner similar to that 
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described in sections 13.21 and 13. 3[ Such a perturbation theory with using a solvabihty 
condition has been employed in many problems [8|, [9l [TO] . 

4. Concluding remarks 

We have formulated a systematic perturbation theory for ([T]). Due to the dilatational 
symmetry ffTIl) . an arbitrary constant A appears in the unperturbed solution (j)u{t)- 
Then, the value of A is determined by the solvability condition fl35l) associated with 
the linear equation ( l33l) for the perturbative correction ip^2^. The advantage of our 
systematic perturbation is in a possibility of developing new and important directions. 
Concretely, following the three problems will be studied soon. 

The first problem is to derive the fluctuation intensity of C{t,t') = 
'^jk '^jit)'^k(t')/NivLst below the mode-coupling transition point for the spherical p-spin 
glass, where aj is a real spin variable that satisfies the spherical constraint Ylf=i '^'j ~ ^■ 
Note that (pit-f) = C{t-t')/C{0) with C(t-t') = (^C{t,t')^ in the stationary regime. 

In a straightforward approach, one may study an effective potential for C{t,t') |llj . 
Indeed, by employing a diagrammatic expansion with neglecting vertex corrections, the 
singular behavior of the effective potential was evaluated [12] . Then, since the minimum 
of the potential corresponds to the solution of ([H), the existence of the dilatational 
symmetry yields the Goldstone mode which carries a divergent part of fluctuations in 
the late stage. More explicitly, A in our expression is treated as a fluctuating quantity, 
and it is identified with the Goldstone mode. (A Related description of fluctuations near 
another type of bifurcation points can be seen in Refs. [131 IE]-) The analysis along 
this line will shed a new light to the understanding of fluctuations near mode-coupling 
transition points. 

As an alternative approach to the description of fluctuations near the mode-coupling 
transition, response properties against an auxiliary external field conjugated to C [15] 
and against a one-body potential field [IE] were investigated. Their studies successfully 
derived the scaling form of a singular part of the fluctuation intensity of C based on an 
idea that such response functions are related to the fluctuation intensity. The extension 
of our work so as to describe the responses may provide a more quantitative result than 
the scaling form. Such an extension is related to a study of dynamical behavior in 
the aging regime, because its behavior is described by a coupled equation of the time 
correlation function C{t,t') and the response function R{t,t'), which are functions of 
two times [H]. In addition to a complicated structure of the equation, the dilatational 
symmetry is replaced with the time reparameterization symmetry. Since the symmetry 
is much wider than the dilatational symmetry, several new features may appear in 
the analysis. See Ref. [18] as a review for the argument on the basis of the time 
reparameterization symmetry. 

The third problem is to analyze a rather wide class of systems with nonlinear 
memory. The qualitative change of the solution f^o of ([2]) is the same type as that 
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observed in an elementary saddle- node bifurcation fT^. Despite this similarity, the 
dynamical behavior near the saddle-node bifurcation is much simpler than that of ([1]) 
owing to the lack of nonlinear memory. Note that an edge deletion process of /c-core 
percolation in a random graph is precisely described by a saddle-node bifurcation [20] 
and it has been pointed out that k-core percolation problems are related to jamming 
transitions [21]. Since nonlinear memory effects might appear in jamming transitions, it 
is important to study a mixed type of dynamical systems which connect the elementary 
saddle-node bifurcation with the mode-coupling transition. The calculation presented 
in this paper may be useful in the analysis of such models. 

By studying these problems, we will have deeper understanding of slow relaxation 
with nonlinear memory. We also hope that this paper provokes mathematical studies 
of the simplest mode-coupling equation ([T]). 
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Appendix A. Approximate expression of cq 

We perform a short time expansion 

oo 

Go{t) = J29nr, (A.l) 

n=0 

which is valid around t = 0. All the coefficients gn can be determined from a recursive 
formula. Concretely, go = 1/2, gi = —1 and g2 = 5/2. Respecting the lowest order 
result Go(t) = go + git, we assume 

= WTW) <^-^' 

for t < t^, where will be determined later. On the other hand, from the asymptotic 
form 

Go{t) = cot-^ (A.3) 

in the limit t — oo, we assume Go(t) = cot^" for t >t^. Since Go(t) is smooth at t = t*, 
we have 

^ = 2co^;^ (A.4) 



l + 2t. 



^ 2acoC-^ (A.5) 



(l + 2t,)2 

These equations lead to = a/[{l — a)2] and cq = (a/(l — a)y{l — a)/2"+^ 
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Appendix B. Derivation of J^2^'^ 

We shall extract a leading order contribution of F2{e^'^'^t2] (pu) in the limit e 0. In 
order to simplify the calculation, we utilize an identity 

r dsfit - s)g'{s) = f ds[f{t - s)g'{s) + g{t - s)f{s)] 

Jo Jt/2 

-f{t)g{0) + f{t/2)g{t/2), (B.l) 

which plays a key role in an efficient numerical integration algorithm for solving mode- 
coupling equations [5]. By substituting (fT6l) into (fT8|) and using this identity, we obtain 

+ g r ds2A'^{s2)[<pl{e-^'\t2-S2)) + 2A^{s2)<p.{t-^'\t2-S2))\. (B.2) 

Jt2/2 

Here, we take At satisfying e''^~"''^ ^ At ^ 1 such that 4>n{^~"''^t2) ^ Ax{t2) in the 
regime At < ^2 < 00. More explicitly, we assume At = e^' with 72 — 71 > 7' > 0. We 
divide the integration regime in the second line of (]B.2[) into two parts, [t2/2,t2 — At] 
and [t2 — At, t2]. Let Ji and I2 be the integration values over the former and the latter 
region, respectively. By a straightforward calculation, we can estimate I2 as 

I2 ~ 2(?A',(t2)(/c + Ax{t2))e^' / dsGo{s)Q{e^'s) 

Jo 

+ g r ds2A'^{s2){Al{t2 - S2) + 2Ax{t2 - ^2)^^(52)) (B.3) 



in the lowest order evaluation. We next combine the second line of flB.3p with Ji and 
return it to the original form. As the result, we obtain 

F.{e-'<H2; 0u) ^ A(t2) - 9{l - fc)Al{t2) +g r ds2Ax\t2 - 52)^/(52) 

^0 

+ 2^7A',(t2)(/e + A,{t2))e^^ / dsGo{s)e{e^^s),iBA) 

Jo 

where higher order terms are ignored. With the aid of ([8]), we rewrite the ffist line of 
(IB. 41) as eA{Xt2) / gc- Furthermore, from an estimation 

dsGo{s)Q{e^'s) ~ coe-^'^^'"^ / dss-^Qis), (B.5) 
Jo 



which is valid in the limit e — > 0, the second line of flB.4p turns out to be of 0(e 



p72-7i(l-a)^ 



These results lead to f l271) . We also find that the higher order terms we have neglected 
in (1B.4P are of 0(e'^2~^^*-^^^"^) by an estimation similar to (]B.5|) . 



Appendix C. Derivation of (1411) 

We take At = e~"', where a' satisfies a' < 71 — 1/2. We also define 

w{ti) = cot^^Q^iti) - Xti. (C.l) 
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Then, for sufficiently small e, h^(t) = 4>e(t) — fc is expressed by 

/i,(t) = Go(t) + 0(e) (C.2) 
for < t < At, and 

K{t) = e^^M^^'t) + e''^[''\e'H) (C.3) 
for At < t < e~'^2 

By substituting 0e(t) = fc + htif) into ([I]), we can write the equation for h^{t). The 

further substitution of flC.2l) and fIC.SP into the obtained equation for yields 

2w'^{t) + 1/8 + 4 / dsi {w{ti - si) - u'(ti)) w'(si) 

= O (^ei/2+7i-«'^ e=^/2^ e^S e"+^/2j . (C.4) 
Extracting the e-independent terms in the limit e — > 0, we obtain 

2w2(ti) + 1/8 + 4 / dsi{w{ti- si) -w{ti))w\si) (C.5) 
We substitute w(ti) = Q{ti) — Ati into this equation. The result becomes (HTI) . 

Appendix D. Short time expansion of Q 

We assume the form 

oo 

Q{h) = Y.qutf'-'^ + \t,. (D.l) 

fc=0 

By substituting (ID.ip into (14T]) . we can determine (A; > 1) recursively from = cq. 
Concretely, the recursion equation becomes 

* = -64c„(V'„' - 1/2) <°-2' 



and 

i , 

Q'fc+i 



2go(M),fc+i - 1/2) 
for > 1, where 



(D.3) 



p(2m + 2n-2) ^ ^ 



with p„ = r(l + an). 
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